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Abstract. Hydro dynamic simulations were performed of the dynamical phase of the 
merging of binary neutron stars (NS-NS) and of neutron star black hole binaries (NS-BH), 
using a physical nuclear equation of state W and taking into account the emission of gravi- 
tational waves and neutrinos. 

Simulations and some results 

NS-NS and NS-BH mergers have been suggested as sources of 7-ray bursts at cos- 
mological distances (e.g., ||) and of r-process nuclei M. These merger events are 
estimated to occur at rates of the order 10 -5 per year per galaxy (e.g., ||). 

By means of three-dimensional hydrodynamic simulations with a Eulerian PPM 
code we have investigated the dynamical phases of the coalescence |l], . The simu- 
lations were performed with four nested cartesian grids which allow for both a good 
resolution near the central black hole and a large computational volume. In a post- 
processing step, we evaluated our models for the energy deposition by vv annihilation 
around the merging objects. 
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Results for some simulations are listed in Table 1. The NS and accretion disk 
(AD) masses are baryonic masses, the spin parameter indicates whether no spin was 
assumed for the NS's ("none"), or the NS spin was taken parallel ("solid") or anti- 
parallel ("anti") to the orbital angular momentum. AA/d is the dynamically ejected 
mass, t s j m the time at the end of the simulation, {L v ) the typical total neutrino 
luminosity, obh = ^ c /(CAfg H ) the relativistic rotation parameter of the black hole, 
q v = L l/ /(Mc 2 ) the conversion efficiency of rest- mass energy of the gas swallowed 
by the black hole into v energy, q uP — E uP /E u the conversion efficiency of v energy 
into e^- pairs, E U v the integral rate of energy deposition by vv annihilation, t acc the 
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Figure 1: Mass- loss phases during NS-NS and NS-BH merging and during the subsequent 
evolution of the massive accretion disk that forms around the BH. 



estimated duration of the accretion by the BH, E u the energy released in neutrinos 
during i acc , and E vj} the corresponding energy converted into e pairs. 

We find that during the dynamical phase of the merging between about 10~ 4 M 
and 10 -2 Mq of cool gas can be ejected (Fig. The exact number depends very 
sensitively on the assumed neutron star spin(s) and on the masses of the merging 
binary components (Table 1). This is of potential interest for heavy-element nucle- 
osynthesis J| . 

Figure pTsketches the result that during the dynamical interaction of the binary 
components (1st phase) low-entropy (dense and cold) matter is ejected from the ends 
of long spiral arms which are formed in the orbital plane. During the merging the 
matter is heated by shocks and viscous dissipation, and intense neutrino fluxes are 
emitted before the gas is accreted into the central black hole. Energy transfer by 
these neutrinos drives a high-entropy gas outflow ( "neutrino wind" ) off the surface of 



short & weak gamma-ray bursts by v v > e~e + 



mass ejection: M ej < few 10 2 M Q 
nucleosynthesis: r-process? 



f ^ M < M 

^\ rem max 



V y nu ^e ar 

NS-NS 
merger jj> 

(~10-»/yr/gal) '"n, ± }f 

M NS = 1.2 ~ 1.8 M e 
M BH = 2~ 10 M 
NS and BH spins 



EOS 



massive NS 
2.2- 2.8 M 
rapid rotation 
v cooling 



BH (Kerr?) 
2 ~ 1 M 

+ 

accretion 
torus 
0.05 - 1 M m 



NS-BH 
merger 

(<10-"/yr/gal) | 




O M NS /M B 



NS mass 
transfer 
cycles 



NO GRB 

UV-burst ? 



At ~ 0.01 ... 1 s 




E T < few 10 51 erg/f n 



NO GRB 

minimum- 
mass NS? 



GW signal: rh < few 10 5 cm 



" ** v signal: < few 10 53 erg 



Figure 2: Possible evolution paths from NS-NS and NS-BH mergers to short 7-ray bursts 
powered by vv annihilation into e pairs. 



the accretion disk. At the same time, energy deposition by vv annihilation in the es- 
sentially baryon-free funnel around the rotation axis powers relativistically expanding 
e ± 7 jets which can give rise to 7-ray bursts. 

Our simulations confirm suggestions that these events could be the sources of 
a subclass of short gamma-ray bursts at cosmological distances (Fig. ||). BH-NS 
mergers yield between 10 and 100 times more energy than NS-NS mergers (Table 1). 
The annihilation of neutrinos and antineutrinos (Fig. |^) which are emitted from the 
hot matter that is accreted into the (forming) black hole, yields total energies E vV 
up to 10 51 erg in our models. It is therefore enough efficient to explain observable 
burst luminosities i 7 ~ E v „ / (fut 7 ) up to several 10 53 ergs _1 for burst durations 
t-y w 0.1-1 s, if the 7 emission is beamed in two moderately focussed jets into a fraction 
fn = 2<5f2/(47r) w 1/100 — 1/10 of the sky. Such a modest amount of beaming (jet 
opening half-angles between about ten and several ten degrees) has to be expected 
because the dense gas in the equatorial plane of the accreting black hole prevents 
relativistic expansion of the pair-plasma jets outside a low-density funnel along the 
system axis. More energy for the burst can be available if Kerr effects of the rotating 
black hole are taken into account. 

The evolution tracks of compact binaries shown in Fig. depend on the NS and 



BH masses, on the total angular momentum in the system, and on the stiffness of 
the nuclear equation of state. Energetic relativistic pair-plasma jets can be produced 
when a hot accretion disk forms around the BH, because high neutrino fluxes are 
emitted over a time much longer than the dynamical timescale of the system. This 
allows for high efficiencies q v and q u p and a large total energy conversion by vv — > 
e + e~. Gamma-ray bursts from these scenarios are accompanied by the ejection of 
nucleosynthesis products in considerable amounts, the emission of large numbers of 
neutrinos with MeV energies, and the production of strong gravitational-wave signals. 
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Figure 3: Azimuthally averaged maps of the rate density of energy deposition by vv annihi- 
lation into e + e~ pairs (in ergcm _3 s _1 , contours represent logarithms, the levels are spaced 
in steps of 0.5 dex) in the surroundings of the accreting BH (white octagonal area at the 
center) which originates from the merging of two NS's with 1.6 Mq baryonic mass (left) or 
of a 1.6 Mq NS with a 5Mq BH (right). In the former case the accretion disk (AD) has 
a mass of 0.26 Mq and the total energy deposition rate below a density of 10 11 gem' 3 is 
4.9 x 10 50 ergs _1 . In the latter case the AD mass is 0.45 Mq and the integral energy de- 
position rate is 15 times higher, 71.6 x 10 50 ergs _1 . The dotted lines represent isodensity 
contours of the AD, the dashed lines show the neutrinospheres for the different v and v 
flavors. The time after the start of the simulation is given in the upper right corner. 



